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We study the interplay between collective and incoherent single-particle motion in a model of two 
chains of particles whose interaction comprises a non-integrable part. In the perturbative regime, 
but for a general form of the interaction, we calculate the spectral density for collective excitations. 
We obtain the remarkable result that it always has a unique semiclassical interpretation. We show 
this by a proper renormalization procedure which allows us to map our system to a Caldeira- 
Leggett-type of model in which the bath is part of the system. 
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Collective motion, i.e., coherent motion of the particles 
in phase space, is a fundamental feature of many-body 
systems. A wealth of information on collectivity is avail- 
able for atomic nuclei [1]. Bose-Einstein condensates 
0] provide more recent examples. There is strong experi- 
mental [F] and theoretical evidence that similar effects 
occur in other fermionic systems as well. Collective mo- 
tioir emerges out of the incoherent, single-particle motion 
whenever favored by energy and kinematic conditions. 
Due to the quantum-classical correspondence principle, 
the collective dynamics on the classical level should be 
reflected in the spectral properties of the corresponding 
quantum mairy-body system. Hence, the spectrum of 
a many-body system comprises states of single-particle 
and of collective character, mixed forms with a partial 
degree of collectivity exist as well. The details strongly 
depend on how the system is probed. 

Consider the Giant Dipole Resonance in heavier nuclei 
as a prominent example which also serves as an inspira- 
tioir for our model to be discussed in the sequel. The 
cross section of electric dipole radiation and the spec- 
tral density of the excitations show at a certain energy a 
huge peak whose spreading width is orders of magnitudes 
larger than the mean level spacing. It can be understood 
in terms of the following picture: the neutrons are con- 
fined to one sphere, the protons to another one. There 
is no or very little relative motion of the nucleons inside 
these spheres. The two spheres, however, move against 
each other, resulting in an enormous response functioir. 
The difference between the center-of-mass coordinates of 
the two spheres is the proper collective coordinate. Many 
other forms of collective motion in nuclei exist. 

Not surprisingly, it is a demanding challenge to under- 
stand the emergence of collective motion and its interplay 
with the incoherent single-particle motion. The vast ma- 
jority of studies in this context relies on effective models 
whose justification is often mainly phenomenological or 
even on the level of hand waviirg if the systeirr iir ques- 
tion is too complex. Better understanding of these issues 
is called for. In the present contribution, we have three 



goals: (i) We want to address the interplay between col- 
lective and single-particle motion froirr first principles in 
the framework of a tractable, yet sufficiently general and 
complex model, (ii) We aim at doing this analytically in 
such a way that we identify the collective coordinate, but 
always keep full control over the single-particle degrees 
of freedom, (iii) We wish to deliver the important insight 
that the spectral density of the collective excitations is 
directly related to classical motion. 

We begin with setting up our model which considerably 
generalizes an integrable model which we studied previ- 
ously 0- In one dimension, two chains a = 1, 2 of in- 
teractiirg particles each with positions x[°'\ i ~ 1, . . . , N 

and momenta p'f'\ i ~ I, . . . , N are coupled to one an- 
other. The total Hamiltonian reads 
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Hq'^^ + Hq''^' is the integrable part 
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model the two chains a = 1,2 before they are coupled. 
The interaction within each chain is harmonic and de- 
scribed by the matrix W which we assume equal in both 
chains. We are mainly interested in selfbound systems 
such as nuclei, where unlike Bose-Einstein condensates 
no external confining potential is needed. It is easy to 
impose corresponding conditions on the matrix W which 
ensure that the interaction is invariant under translations 
and the system is bounded 0. We now couple the two 
chains by an interaction which depends on the differences 
between their coordinates. 
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For every choice of the coupling matrix K, the Hamilto- 
nian Hq is translation invariant. Clearly, the model is up 
to now integrable. 



2 



We generalize the model by adding the translation in- 
variant term A-ffi which breaks integrability, 



N 



,(1) 
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where / is an arbitrary, positive and even analytical func- 
tion of the form 



f{z) = J2 /"^ 



(5) 



Tl = 2 



We introduce the parameter A, because we aim at a per- 
turbative discussion. 

To quantize our model, we replace coordinates and 
momenta by operators i^"'' and Importantly, we 

make this step on the level of the original particle de- 
grees of freedom. Motivated by the above mentioned 
Giant Dipole Resonance, we aim at studying collective 
excitations and the associated spectral density. We ex- 
pect that it shows a pronounced peak, which we wish to 
understand in (semi) classical terms. Naturally, the col- 
lective coordinate X is the difference between the mass 
centers of the two chains, and it is convenient to rescale 
it with a factor N/2, 
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and accordingly for the collective operator X. To probe 
the existence of quantum collective states in the excita- 
tion spectrum, we investigate the correlator 



S{t) = ($o|^(t)^(0)|$o) 
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where |<I>o) is the ground state of the total Hamiltonian 
H. Here, X{t) is the Heisenberg picture of the operator 
X with the time evolution governed by the total Hamil- 
tonian H. The Fourier transform of the correlator ([7]), 



ti=0 



Efj, — Eq 
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is the desired spectral density of the collective excita- 
tions. It measures the strength of the transitions between 
the ground and the excited states \^fj} of the whole sys- 
tem and can be interpreted as the response of the system 
that is excited by the transition operator X. Following 
the terminology in many-body physics, we say that there 
is a collective quantum state for an energy Ecoi = Eq + Huj 
if S{uj) (smoothened over some energy interval) has a pro- 
nounced spike at the corresponding frequency uj. 

To leading order A in perturbation theory, we obtain 



the following expression for the correlator 
S{t) ^ {0\X I {t)X I {0)\0) 
+ AE {al{0\Xjit)Xjmi) + ialr{l\Xjit)XjmO)) 
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dh{0\[Hi{ti),Xi{t)]Xi{Q)\0) + OiX") , (9) 



where the sum runs over the eigenstates |/) of Hq, and 
where |0) is the ground state of Hq. The coefficients ajg = 
{l\Hi |0) / (Eq — Ei) turn out to be real due to time reversal 
invariance. The sum in Eq. ([9]) arises from the correction 
to the ground state while the last term results from the 
perturbation of the Hamiltonian. The collective operator 
Xj and the non-integrable part Hu of the Hamiltonian 
appear in the interaction picture whose time evolution 
is governed by the integrable Hamiltonian Hq. For any 
operator F, we have 



Fi{t)^exp[ -Hot]Fexp 
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For later purposes, it is useful to consider the imaginary 
part of the correlator, ImS'(t) = Si{t). We notice that, 
by virtue of Eq. ([S]), the Fourier transforms of S{t) and 
Si{t) are connected through 



S{lu) = 2ie{uj)Si{Lu) 



(11) 



where Q(uj) denotes the Heaviside step function. Since 
lm{0\Xj{t)Xj{0)\l) vanishes for I ^ 0, the imaginary part 
of the correlator simplifies, and we find 

Si{t)^lm{0\Xiit)XiiO)\0) 



dh{0\ [Hu{ti),Xiit)],Xi{0) \0)+O{\'') 



(12) 

At this point, we may use our previous results The 
first term of Eq. (|12p was evaluated by mapping Hq into 
the form of a Caldeira-Leggett-like model, in which X 
can be viewed as the coordinate of a "big" particle in 
a harmonic potential which is coupled to a "bath" of 
harmonic oscillators. The interpretation of the "bath", 
however, differs significantly form the standard Caldeira- 
Leggett situation [8|. In our case, the "bath" is not ex- 
ternal, it is part of the system and formed by the inter- 
nal degrees of freedom. The resulting expression for the 
spectral function S{io) is then 
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— Lu^' — iuj^iuj) 

where 7 formally coincides with the classical "damping" 
kernel, but here it describes the spreading of the collec- 
tive excitation over the spectrum. There is not an en- 
ergy loss or any kind of dissipation in our system. The 
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resonance frequency f^o is the fundamental oscillator fre- 
quency of the corresponding classical problem. To il- 
lustrate this result, we mention that, when the collec- 
tive degree of freedom X interacts with an ohmic "bath" 
(see Ref. the spreading kernel 7(w) = 70 is a con- 
stant and S{uj) has a Lorentzian shape with the width 
7o at the position of the classical oscillator frequency 

no = Vf^o - (70/2)2. 

We now consider the crucial second term on the right 
hand side of Eq. (fT2)) . A naive continuation of the ap- 
proach in Ref. [Tj] quickly becomes cumbersome. Luckily, 
there is much better way of tackling this term which even- 
tually leads to a new insight to our problem. The second 
term is easily seen to be of the form ^ dtix{ti, t), with 
the kernel x(ii,i) given by 
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Owing to the harmonic form of Hq, we find the identity 
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for the commutator of the n-th powers of differences with 
the collective operator. For all n, it can be reduced to the 
operator (a;-^'' — x'"j^^)"'~^ multiplying the function l3ij(t) 

'xf^),X{-t)] = I3ij{t)l. We notice that 



defined by [(: 

the commutator in the latter expression is proportional 
to the unit operator 1. Applying this formula twice yields 



(16) 



where the elements of the matrix C are given by the 
ground state expectation values involving the second 
derivative of the arbitrary function / defining the non- 
integrable part of the interaction in Eq. ^ , 
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We emphasize that this result is not due to an expansion 
of the function f{z), it applies in leading order A to all 
functions of the form ([5]). 

We arrive at the important insight anticipated above: 
precisely the same equation for the kernel xi^ii t) follows 
when using the harmonic Hamiltonian 
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The Cij themselves depend on the parameters of the 
integrable Hamiltonian Hq. Starting from the definition 
([T7|. we derive 



Cij — 



+00 

^//"(xyZi^)exp(-;^ 
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xf^f\0) can be 
They are given by 




where the quantities a.y — {0\{x 
related to the parameters of Hq 
aij — (Tii + Tjj)/A, where is the fc-th diagonal ele- 
ment of the matrix 



r = 



and Af is a diagonal matrix whose elements read 
Mij — Sij J2i ^u- The coefficients Cij depend on h since 
they result from sandwiching the function / with the 
ground states of the integrable Hamiltonian. If we re- 
strict ourselves to leading order ft, only the first term of 
the expansion ([SJ enters. In such a semiclassical regime 
Eq. simplifies and we have Cij = 8/2 (F^i -I- Tjj). 

Since H^ is harmonic, the spectral function S'(w) can 
now be calculated using Eq. (I13|) . The renormalized oscil- 
lator frequency Vl^ and the renormalized spreading ker- 
nel 7^ are determined by the renormalized coupling con- 



Kij -\- \Cij/2 and by Wij. The spread- 



stants K^, 

ing kernel can be expressed through the spectral density 
function (T^ (a;) [9.], 
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Employing our previous results Q, we express the spec- 
tral density function through the parameters of the har- 
monic Hamiltonian. This yields 



a«(c.) = - 



-Imk 



2TTmuj {uj + ie)l - [2Krlm)^/' 
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where Kr and k are obtained from the matrix = 
A^{W + -f K^)A, with A being a discrete N x N 
cosine transform (DCT), see Ref. [7|. More precisely, 
Kr is the (A^ — 1) X (TV — 1) matrix obtained from 
by deleting the first row and the first column while the 
vector k is the first column (excluding the first element) 
of K^. The renormalized oscillator frequency turns out 
to be 



In other words, the effect of a general, non-integrable 
perturbation can, to leading order A, be fully accounted 
for by a proper renormalization of the integrable Hamil- 
tonian Hq. Since H^ is harmonic, the spectral function 
S{u)) is given by Eq. (IT51) . where the renormalized oscil- 
lation frequency and the spreading kernel 7^ depend 
on the constant matrix elements Ca. 



(23) 



It is important to notice that the spreading kernel (|2ip . 
the oscillator frequency ([23]) and therefore the spectral 
density S{uj) are fully determined by the classical dy- 
namics of the renormalized Hamiltonian H^ . In partic- 
ular, using our previous results J}, the equation for the 
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classical time evolution of the collective coordinate reads 
^ + in^mt) + j ,^it - s)^ds ^ . (24) 



The same equation holds for the expectation value {X) (t) 
if the initial state of the system is properly chosen. Equa- 
tion ([M)) describes a damped, i.e. in our case spread, har- 
monic oscillator whose spreading kernel 'y^{t) and oscilla- 
tor frequency Hq are given by Eqs. (PT|) and (1^ . respec- 
tively. Importantly, the solution of Eq. ([^ also deter- 
mines, for properly chosen initial conditions, the spectral 
density of the collective excitations. We notice, however, 
that itself contains quantum corrections which de- 
pend on h. Put differently, is identified as the proper 
effective Hamiltonian whose classical dynamics — rather 
than the classical dynamics of the original, total Hamil- 
tonian H — determines the spectrum of collective exci- 
tations in leading order A. We also notice that higher 
order terms in the perturbative treatment of S{t) come 
with higher powers of h. Indeed, it is straightforward to 
see that in the case of f{z) = f2z'^ the n-th term of the 
perturbative expansion scales as (ftA)". In this sense the 
renormalized Hamiltonian provides the first semiclas- 
sical correction to the spectrum of the collective modes. 

We briefly discuss the conditions for the validity of 
our perturbative approach. The approximation ^ can 
be used if the following conditions are satisfied, see 
Ref. [13] • First, the gap between the ground state and 
the first excited state of Hq must be sufficiently large, i.e., 
A|(0|i/i|0)| ^ Swmin, where Lorain is the minimal oscilla- 
tor frequency of the classical system given by the lowest 
eigenvalue of the matrix W + M + K. Second, the time t 
of propagation must be bounded by X\{0\Hi\0)\t/h <^ 1. 
As we are interested in time scales of order t ~ ^q^^ 
the first condition implies the second one. Under these 
conditions the spectral characteristics such as energy and 
spreading width of the collective excitations are close to 
their values for the unperturbed Hamiltonian Hq. We 
notice, however, that such a "small" perturbation XHi 
might be quite large on the scale of the mean level spac- 
ing which is of the order . This means that the local 
distribution of the energy levels of the Hamiltonian H 
might be essentially different from the energy level dis- 
tribution of the integrable Hamiltonian Hq . 

In conclusion, we studied, in the framework of a sim- 
ple model, the emergence of collectivity from first princi- 
ples. We did not start from an effective model, we rather 
derived an effective description and still kept full con- 
trol over the original degrees of freedom. In doing so, 
we relate the expectation value of the collective opera- 
tor and the spectral density of the collective excitations 
to a purely classical equation. We consider that to be 
important, as it can be viewed as a justification of the 
routinely used strategy in many-body physics, where ef- 



fective models are set up classically and then quantized. 

Beside the fundamental aspects just mentioned, there 
is further motivation for our study: Statistical analysis 
of the spectra indicates that collective motion is typ- 
ically regular while the incoherent single-particle mo- 
tion yields spectral statistics described by random ma- 
trices, see Refs. ill- 13 1 . This coexistence of both regular 
and chaotic dynamics in the same system is a truly in- 
triguing dynamical aspect of many-body systems [l3 |. 
The regularity of collective motion implies that the re- 
cent arguments [151] strongly supporting the Bohigas- 
Giannoni-Schmit conjecture for single-particle systems 
do not carry over in a straightforward manner to many- 
body systems. In short, this conjecture states that 
the spectral statistics of the single-particle system is of 
random-matrix type if the corresponding classical system 
is chaotic. Our study is thus needed when addressing the 
role of collectivity in quantum chaos. 
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